Consequences of energy conservation in relativistic heavy-ion collisions 
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Complete characterization of particle production and emission in relativistic heavy-ion collisions 
is in general not feasible experimentally. This work demonstrates, however, that the availability of 
essentially complete pseudorapidity distributions for charged particles allows for a reliable estimate 
of the average transverse momenta and energy of emitted particles by requiring energy conservation 
in the process. The results of such an analysis for Au+Au collisions at ^/sMN= 130 and 200 GeV are 
compared with measurements of mean-pT and mean-iSr in regions where such measurements are 
available. The mean-pT dependence on pseudorapidity for Au+Au collisions at 130 and 200 GeV is 
given for different collision centralities. 

PACS numbers: 25.75.-q, 13.85.Ni, 21.65.+f 



INTRODUCTION 



In recent studies of ultra-relativistic heavy-ion colli- 
sions at the Relativistic Heavy-Ion Collider (RHIC) at 
Brookhaven National Laboratory much attention has 
been devoted to measurements of particles and radia- 
tion emitted in the transverse direction. This focus is 
well justified because particles emitted in this direction 
are expected 0] to carry information about the region 
of highest energy density formed in the process, which 
is believed to consist of de-confined quarks and glu- 
ons M- Experimental evidence based on Au-|-Au colli- 
sions |3, 0, IE 01 indicates that a phase of matter with ex- 
treme energy density is indeed formed in such collisions, 
the properties of which correspond to an ideal fluid of 
quarks and gluons. 

The present paper is instead focussed on following the 
total energy in the collision and how it is distributed in 
the exit channel. It is interesting to note that the number 
of particles emitted at mid-rapidity increases only slowly 
as a function of collision energy l7| ■ The increase in col- 
lision energy by more than a factor of one hundred from 
the CERN SPS facility to RHIC thus result in only dou- 
bling of the mid-rapidity particle multiplicity. The bulk 
of the additional energy is carried off by particles emit- 
ted away from this region. Thus, in order to study the 
energy balance in the reaction it is necessary to measure 
particles emitted over essentially the full solid angle. The 
PHOBOS experiment at RHIC 8] is unique in being close 
to achieve such coverage for charged particles and nearly 
complete distributions of charged particles iiipseudora- 
pidity space, dNch/drj, have been reported In this 
work, these and other data are used to show that energy 
conservation places severe constraints on the transverse 
momentum distributions away from midrapidity, outside 
the range where they can be measured. The results are 
used to calculate the pseudorapidity dependence of the 
mean transverse momentum, mean transverse energy and 
an estimate of the initial energy density in the spatial 
region from which particles with large longitudinal mo- 
menta originate. 



II. ENTRANCE CHANNEL 



At relativistic energies, it is normally assumed that 
only a fraction of nucleons, Npart, of the incoming heavy 
ions take part in the collision. The non-participaring 
nucleons, the spectators, are assumed to continue un- 
affected by the collision between the participants and 
do therefore not contribute to the particle production. 
For the purpose of the present work, the energy of the 
spectators will therefore not be included in the energy 
accounting. At the 130 GeV and 200 GeV collision ener- 
gies studied here, it is believed that the pseudorapidity 
distribution of the spectators falls outside the PHOBOS 
acceptance. On the other hand, it is assumed that all of 
the energy carried by the participants will appear in the 
exit channel. In a recent analysis of net baryons in the 
rapidity range y=0-3.2, the BRAHMS collaboration 
have concluded that they retain ~27% of their initial en- 
ergy, which might lead one to the conclusion that only 
~ 73% should appear in the exit channel. It is impor- 
tant, however, to realize that also the partially stopped 
net baryons are included in the measured (and extrapo- 
lated) dN/drj distribution of charged particles considered 
in this work as demonstrated in Sect. VII. Consequently, 
one needs to account for all of the energy of the partici- 
pants, namely 



^/SNN X Npart/2, 



(1) 



where ^^snn is the center-of-mass energy per nucleon 
pair and Npart/ '2 is the number of such pairs. The value 
of Npart is normally calculated in the Glauber model 
for collisions at different impact parameters using exper- 
imental nucleon-nucleon cross sections. The centrality of 
heavy-ion collisions is typically given by the fraction of 
the total cross sections such that the 0-6% centrality bin 
corresponds to essentially head-on collisions exhausting 
6% of the total cross section, the 6-15% centrality bin is 
associated with the 9% of the cross section with increas- 
ing impact parameter, etc. 
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III. EXIT CHANNEL 

Since the analysis is based on measurements of pseu- 
dorapidity distributions of charged particles, dNch/dr/, 
one should expect to account for only this fraction of the 
total energy, i.e. 

= fch X X Npart/2, (2) 

where fch is the fraction of exit particles that are charged. 

The energy of a particle is E = ttit cosh?/, where ttit = 
\J TV? + is the transverse mass, m is the rest mass, pt 
is the transverse momentum, and y is the rapidity of the 
particle. Using the identity mrsinhy = pTSvahrj one 
finds 



E = rriT cosh.y = y inn? + p^cosh^ "q. (3) 

The total energy of charged particles may be evaluated 
by summing over pions, kaons, and protons and integrat- 
ing over the distribution in px and 77 

where 

dE ^ r n — I — d^^^ , 

— = > Pi Wmf+p^cosh 77- — —dpT, 5 
dr] -r^ Jo ^ dprdr] 

and Pi is the relative abundance of particles with mass 
rrii. 

In the Appendix, the effects of computing dE/drj (Eq. 
5) under different assumptions are explored and it is 
shown that using an average mass (m) — J2i=Tr k p Pirni 
is an excellent approximation to an accurate evaluation 
of the individual contributions from the three particle 
species with different (p^) values. The dependence of 
(pt) on particle mass has been observed experimentally 
and is normally interpreted as a consequence of radial 
flow of the hot fireball. In the following, the expression 




is therefore used to find those ry-dependences of {pr) that 
satisfy energy conservation. 



IV. INPUT DATA 

The dNch / drj distributions the and corresponding val- 
ues of Npart were taken from Ref. . The data and error 
bands were fitted to a high accuracy using three Gaus- 
sians, one centered at 77 = 0, and two of identical shape 
centered at ±770 using five fit parameters. 



A. 130 GeV 

Measurements by STAR jH El Ei have shown that 
the overall charged particle pr-distributions near mid- 
rapidity follow a power-law behavior, 

dpT Pt \ PoJ 

where n and Pq are adjustable parameters. The values of 
n and (pt) = '^Po/i^ — 3) were taken from |l4l | and in- 
terpolated to match the centrality bins of the PHOBOS 
dN/drj data. In the Appendix it is shown that the total 
energy in the exit channel is very insensitive to the par- 
ticle rest mass. This power-law shape of the px spectra 
(in 77-space) is therefore used for all three particle types 
(pions, kaons, protons/anti- protons). 

The charged particle fraction was estimated using iden- 
tified particle multiplicities dN/dy\y^o from Rcfs. E3| 
(pions), E3 (kaons), and 16] (protons/anti- protons) and 
converting to pseudorapidity space dN/dr]\ri=o by using 
the fitted spectral shapes. The results were again inter- 
polated to match the PHOBOS centrality bins. 



B. 200 GeV 

As yet, no analysis of charged hadron py-spectra has 
been published for 200 GeV Au-|-Au collisions. In this 
case it is therefore assumed that the centrality depen- 
dence of the 7^ parameter is identical to that found at 130 
GeV. In order to estimate {pT)ri=o and fch for uniden- 
tified hadrons at this energy, the values of dN/dy\y=o 
of charged particles were taken from STAR Uj and 
PHENIX 18], converted to pseudorapidity space using 
the fitted spectral shapes, interpolated to match PHO- 
BOS centrality bins and averaged (with equal weights). 
Note that the transformation from rapidity to pseudora- 
pidity space increases the values of (pt)))=o to be larger 
than those obtained directly from the dN / dpT-spectTa 
measured for a fixed rapidity interval. 

The charged particle fractions, fch, were obtained from 
the same data according to the method used for 130 GeV. 
All resulting values are listed in Table I. 

V. AVERAGE TRANSVERSE MOMENTUM 

In this section it is shown that the requirement of en- 
ergy conservation and the fact that (px) is experimen- 
tally known at 7/ = imposes severe restrictions on the 
PT-distribution of particles emitted at large pseudora- 
pidities. 

First, it is shown that a constant (as a function of 77) 
(pt) does not conserve energy. By assuming that par- 
ticles emitted at all rj values are characterized by n= 
21.9 and (pt) = 0.544 (GeV/c), corresponding to those 
measured at mid-rapidity for 200 GeV Au-f Au collisions 
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FIG. 1: Panel a) Pseudorapidity distrbution for 0-6% central 
Au+Au collisions at 200 GeV |3| (solid points and error band). 
Panel b): dE/dr) is plotted for different assumptions of the 
77 dependence on (pt), which are shown in Panel c). The 
estimates of the r;-distribution of net-protons and their energy 
derived in Sect. VII is shown a grey shading in panels b) and 
c), respectively. 



at 0-6% centrality, Eq. 6 gives a total energy of E^f^— 
35,523 GeV, which is substantially larger than the £^°^* 
= 21,844 GeV expected in outgoing charged particles on 
the basis of energy conservation. One must therefore 
conclude that (pt)('7) assumes smaller values (softer px- 
spectra) away from mid-rapidity. In the following, dif- 
ferent shapes of the (pT)-dependence on 77 are explored. 
Guided by predictions of theoretical models, such as HI- 
JING Hi, RQMD [13, and AMPT 13 it is assumed 
that {pt) falls off monotonically away from 77=0. Only 
analytical functions which have this property are there- 
fore considered. 

Assuming a Gaussian 77 dependence of (pr) and requir- 
ing energy conservation allows for a determination of the 
standard deviation a of the Gaussian dependence. The 
solid curves in Fig. 1 represent a calculation of dE/drj 
for 200 GeV central (0-6%) Au+Au assuming 



(pT)oexp(-r7^/2cr^ 



(8) 



using the parameters {pt)o — 0.544 GeV/c, cr=4.53. 
These parameters result in a total charged particle energy 
which reproduces the expected value of E°J^* = 21,844 
GeV (obtained from Eq. 2 using parameters listed in 
Table I). 

Although a Gaussian ?7-dependence of the mean trans- 



verse momentum probably approximates the actual de- 
pendence quite well other dependencies have also been 
studied. One parametrization is given by 



1 + exp(-?7g/a^) 
° 1 + exp[(r7 - 770) (77 + ?7o)/a2 



(9) 



i.e. SL flat-top distribution which falls off at ±770 over a 
range of a. Since this form contains two adjustable pa- 
rameters, a family of curves exists that will obey energy 
conservation as illustrated in Fig. 1 where such curves 
are shown for a = 2 (dashed) and a = 3 (dotted-dashed) 
for the 200 GeV 0-6% central Au+Au coUision data. 
Finally, a triangular" shape given by 



{pt)t, = {pt)o 



1 - |?7|/a- bexp(-|?7|/b)/a; 
l-b/a 



(10) 



where a controls the slope away from 77 = and b the 
radius of curvature at midrapidity is explored. This func- 
tion does not deviate significantly from the Gaussian 
when both are required to yield the expected total en- 
ergy. 

130 GeV Au+Au 
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FIG. 2: Mean transverse momentum is shown as a function of 
pseudorapidity for eleven different centrality bins in Au+Au 
collisions at 130 GeV. Solid curves are Gaussian pseudora- 
pidity depencences whereas dotted curves represent flattop 
dependencies, which also conserve energy during the collision. 
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It is interesting to note that the four curves shown in 
Fig. Ic representing different 77 dependences of (pt), all 
converge at at a value of (pt) = 0.312 GeV/c at 77 = 
5.05. This value of (pr) value may be considered a direct 
experimental measurement because it is independent of 
the assumed functional form of the (pt) -dependence on 
■q. 




on the charged-particle fraction, Sfch, which has been 
propagated from the measured values of dN/ dy for identi- 
fied pT spectra of pions, kaons and protons and converted 
to the pseudorapidity space using the relevant pt fit func- 
tions. The second term, SNpart, was taken directly from 
Ref. 0, as was the last component SNch- The final er- 
ror on the (pt) {rj) functions were computed as a function 
of pseudorapidity. Clearly, the contribution from S{pt) 
dominates near 77 = 0, whereas the errors on dNch/drj 
and prevail at larger values of jTyj. In comparison to 
the uncertainty on experimental input data, the errors 
stemming from using an average rest mass of all hadrons 
(too ~ 0.208 GeV/c^) are insignificant as demonstrated 
in the Appendix. 

200 GeV Au+Au 




FIG. 3: Same as Fig. 2, but for 200 GeV. 



FIG. 4: Comparison of mean-pT obtained from different 
methods. See text for details. 



VI. RESULTS 

The Gaussian (pT){r/) functions, consistent with en- 
ergy conservation, arc shown for 130 and 200 GeV 
Au-f Au collisions in Figs. 2 and 3 for six different cen- 
trality bins in each case. The parameters describing these 
functions are listed in Table I. Based on pT spectra, 
measured by BRAHMS [HIS^ over a limited (pseudo)- 
rapidity range it is believed that the Gaussian shape is 
realistic. 

The error bands were obtained by propagating the er- 
rors on the experimental input parameters used in the 
analysis. These are: d{pT),SEch = \/{5fchY + {SNpart}'^ 
and SNch- The uncertainty on the mean pt of hadrons, 
Sipr), was obtained from Refs. [11 [ll[l3. The un- 
certainty on the total charged particle energy in the exit 
channel, SEch, has two sources, namely the uncertainty 



VII. COMPARISON TO BRAHMS RESULTS 

In recent publications, the BRAHMS collaboration 
has presented transverse momentum spectra and ra- 
pidity densities for pions, kaons |22l | and protons/anti- 
protonsfToj for rapidities ranging from y — to y ^ 3.4. 
It is therefore of interest to compare these data to the re- 
sults of the present analysis. Using double-gaussian fits 
to the BRAHMS data on (pt) and dN/dy, and by using 
the spectral functions (exponential for pions and kaons 
and Gaussian for protons/anti-protons) for the shape of 
the pT-spectra, the values of (pt)(??) has been calcu- 
lated over the full pseudo-rapidity range. The results 
are shown as solid circles in Fig. 4 in the range where 
experimental data are available and as a solid curve in 
the extrapolated region. The BRAHMS data fall within 
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TABLE I; Summary of relevant parameters used in the analysis are listed as a function of the collision centrality. The origin 
of the input data (Npart, fch,n, ajid{pT)n=o) is detailed in the text. The last column list the a of the Gaussian {pT){'n) curve 
(Eq. IHl that are consistent with the value of Etot- 



y/s]^=130 GeV 


Gaussian 


Bin 


part 


fch 


n 


TT'OUt 

Ech 




a 










(GeV) 


(GeV/c) 




0-6% 


340±11 


0.627±0.017 


22.1 


13857±585 


0.519±0.008 


4.34 


6-15% 


275±9 


0.627±0.017 


20.7 


11208±476 


0.516±0.008 


4.04 


15-25% 


196±8 


0.627±0.020 


18.5 


7988±414 


0.511±0.008 


3.68 


25-35% 


136±6 


0.625±0.023 


17.2 


5525±317 


0.504±0.007 


3.42 


35-44% 


90±5 


0.624±0.029 


16.4 


3650±264 


0.490±0.008 


3.32 


45-55% 


59±5 


0.622±0.034 


14.8 


2385±208 


0.478±0.007 


3.36 


^31^=200 GeV 


Gaussian 


Bin 


part 


fch 


n 


TT'OUt 


{Pt)v=o 


a 










(GeV) 


(GeV/c) 




0-6% 


344±11 


0.635±0.025 


22.1 


21844±1108 


0.544±0.021 


4.53 


6-15% 


274±9 


0.636±0.026 


20.7 


17426±912 


0.538±0.021 


4.42 


15-25% 


200±8 


0.636±0.027 


18.5 


12720±742 


0.530±0.021 


4.21 


25-35% 


138±6 


0.638±0.029 


17.2 


8804±554 


0.519±0.021 


3.82 


35-45% 


93±5 


0.639±0.030 


16.4 


5943±424 


0.508±0.021 


3.78 


45-55% 


59±5 


0.640±0.031 


14.8 


3776±315 


0.496±0.021 


4.06 



the error band of the present analysis, but lie in the up- 
per range near mid-rapidity. The open circles represent a 
triangular (a=11.2, b=1.2) 77-dependence, which exceeds 
energy conservation by 13%, but follows the extrapolated 
BRAHMS data quite closely. The open triangles achieve 
energy conservation and are seen to follow the Gaussian 
shape obtained in the analysis presented in the previ- 
ous section (note that the slightly higher value at 77=0 
accounts for an insignificant portion of the total energy 
balance). 

In a recent analysis, the BRAHMS collaboration has 
found that ~ 27 ± 5% of the initial energy remains in 
net baryons after the collision]^. In agreement with 
predictions by the HIJING/BB model Q it is found 
that the net-protons are emitted in the region ?/ ~ 4± 2, 
which explains why a relatively small number of parti- 
cles can account for such a large fraction of the available 
energy. A simple model which approximates the rapidity 
distribution of net-baryon predicted by Ref. as 



dN 
dy 



— = 2.2coshy[e"^™''^"^^^ 6"^'''"''^+^^^ Z^^*^"] 



(11) 



and uses an exponential p^— distribution, dN/dpr — 
pj,e-PT/t ^iti^ ^ ^ 0.5e-f'/'^2 (MeV/c) also accounts 
well for the predicted rapidity distribution of net- 
baryon energy per participant, dE/dy, see Fig. 5. This 
model also predicts the pseudorapidity distribution of net 
baryons and their energy shown as dashed curves in Fig. 
5. The pseudo-rapidity distributions of net protons and 
the energy they carry according to this simple model are 
also shown as shaded regions in Figs, la and lb, re- 
spectively. It is evident that the predicted net proton 
distribution is easily accomodated within the measured 



dNch/drj distribution. The associated energy appears, 
however, to exclude a range of (pt) dependencies using 
the flattop fimction with a < 3. Functions with a sharper 
maximum at 77 = 0, such as the Gaussian or triangular 
functions studied here, are clearly preferred. 
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FIG. 5: Panel a): rapidity distributions of net baryons are 
shown for HIJING/BB (solid histogram) and an analytical fit 
using Eq. 11 (solid curve). The dashed curve represents the 
pseudorapidity distribution corresponding to the fit. Panel 
b): The corresponding energy distributions are shown. 
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VIII. COMPARISON TO THEORY 

The resuhs concerning the {pr) dependence on pseu- 
dorapidity obtained in the previous section were based 
almost exclusively on constraints of experimental data; 
the measured dNch/drj and net proton distributions in 
pseudo-rapidity space. It is, however, of interest to com- 
pare these results to predictions of theoretical models of 
particle production in relativistic heavy-ion collisions. In 
a recent work, Topor-Pop et al. !24'| have presented the 
predictions for the transverse energy of the HIJING/BB 
and RQMD |0| models. 



1 30 GeV Au+Au 0-6% 



400 




-2 2 
Pseudorapidity, r| 



FIG. 6: Comparison of charged-particle transverse energy 
with predictions from the HIJING/BB model. See text for 
details. 

The mean transverse energy of charged particles as 
a function of pseudorapidity obtained in the present 
analysis using a Gaussian 77-dependence of (pr) (dashed 
curve) is in Fig. 6a compared to those predicted by 
the HIJING/BB (thick histogram) and RQMD (thin his- 
togram) models. The overall distribution of transverse 
energy obtained in the present analysis agrees reason- 
ably well with the model predictions. The value at mid- 
rapidity exceeds the prediction somewhat, which is also 
in agreement with the measurement of Ref. j2fii |. The 
latter has been multiplied by a factor of /cft,=0.627 to 
account for the fact that only charged particles are con- 
sidered here; the PHENIX experiment refers to both 
charged and neutral particles. Fig. 6b shows a compari- 
son of the charged particle transverse energy per charged 
particle, i.e. {dE^ / drf) / [dN ch / drf) between the present 



analysis and the HIJING/BB and RQMD predictions. 
While the RQMD model predicts a rather flat depen- 
dence on pseudorapidity, we find that the nearly trian- 
gular shape derived from the present analysis is quite well 
reproduced in the HIJING/BB calculations. The closed 
and open circles in Fig. 6b represent values obtained us- 
ing the PHENIX j26] and PHOBOS data on dNch/drj, 
respectively. The latter is in better agreement with the 
present analysis (dashed curve), which is based on the 
experimentally derived value of (pt)?7=o listed in Table I. 



IX. ENERGY DENSITY 

A crucial prerequisite for achieving the quark-gluon 
plasma phase is the initial energy density produced in 
the collision. Based on lattice QCD calculations [23| it 
is estimated that the transition to this phase occurs at 
an energy density of etr ~ 0.7-1.0 (GeV/fm'^). Differ- 
ent estimates of the maximum energy density reached in 
central Au-I- Au collisions at RHIC energies 0, all in- 
dicate that the predicted threshold value of etr in these 
collisions is exceeded by a large factor dependent on the 
assumed value of the equilibration time tq, for which a 
conservative value of tq = I (fm/c) is often used. Esti- 
mates of the energy density are normally based on the 
Bjorken prescription jlj] applied to particles emitted near 
mid-rapidity that probe the highest energy density of 
the fireball. It is, however, also of interest to use the 
same prescription to study how the energy density falls 
off away from the central region. Since the longitudinal 
position z of a particle at the time tq is related to its 
rapidity y by z — tq sinh y one can study the initial en- 
ergy density felt by particles as a function of rapidity, 
and therefore, pseudorapidity. Starting from the normal 
Bjorken estimate 



rriT dN 
TqA dy ' 



(12) 



where A is the overlap area of the two nuclei, one may 
write 



€0(77) 



I dN 
TqA drj 



(13) 



i—7r,K,p 



mf + p\ cosh 77 fif^ 



Pt cosh 77 



dpi 



dpi 



by transforming into 77-space. 

The resuh is illustrated in Fig. 7 for 200 GeV Au-|-Au 
collisions at 0-6% (solid curve) and 0-40% (dashed curve) 
centrality. The quantity cqTo is seen to be sharply peaked 
at mid-rapidity and falling off rapidly for larger values of 
\ri\. Note that this result is about 25% larger than ob- 
tained in other works. The reason is that the overlap 
area. A, for 0-6% is actually ~25 % smaller than for a 
head-on collision, which has often been used. As a simi- 
lar behavior has been seen in the pseudo-rapidity depen- 
dence of the elliptical flow signal, V2, measured for 



7 



"1 ' I ' I ' r 
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FIG. 7: Energy density, eo times the equilibration time, tq cal- 
culated for 200 GeV Au+Au collisions at 0-6% (solid curve) 
and 0-40% (dashed curve) centrality, the latter corresponding 
to the centrality range of the elliptical flow measurements of 
Ref. [2^ . is shown as a function of pseudorapidity. A Gaus- 
sian (pt) dependence on rj was used which satisfies energy 
conservation. 



0-40% centrality. It is reasonable to assume that the ideal 
liquid behavior observed at mid-rapidity via the strong 
flow signal is associated with the high energy density in 
this region. Consequently, it is tempting to speculate 
that there is a direct relation between the 77-dependence 
of V2 and that of the energy density eg. In this context 
it is noteworthy that purely hydrodynamical models p9| 
have had difficulty reproducing the observed pseudora- 
pidity dependence of V2, although recent work [s^l ap- 
pears to have been more successful. 



X. CONCLUSIONS 

By requiring energy conservation in ultra-relativistic 
heavy-ion collisions it has been shown that unique in- 
formation about the average transverse momentum of 
charged particles emitted in such collisions can be ob- 
tained in regions of phase space, where direct experi- 
mental measurements are not available, and, indeed, are 
virtually impossible. Based on (px) measurements over a 
limited rapidity region y = to y = 3.4 by the BRAHMS 
collaboration, it is believed that a Gaussian r] dependence 
is quite accurate. Under this assumption, predictions for 
(pt) are given for Au + Au collisions at ^sjvjv=130 and 
200 GeV as a function of centrality. In a recent analy- 
sis of the net-proton abundance away from mid-rapidity 
by the BRAHMS collaboration ^ it is concluded that 
only about 73% of the available collision energy leads to 
particle production, with the remaining energy being car- 
ried away by net baryons leaving the collision with near 
beam rapidity. Based on predictions of the HIJING/BB 



model the present analysis shows that these particles are 
included in the measured (iiVcft/c''7-distributions. Thus, 
the full energy of the participating nucleons is needed 
to account for the observed charged-particle multiplicity. 
Based on the Bjorken prescription, the pseudorapidity 
dependence of the energy density has been derived as- 
suming an ?7-independent equilibration time, tq. For 200 
GeV central (0-6%) Au-|-Au collisions the mid-rapidity 
energy density peaks at ~ 6.5/to (GeV/fm^), but falls 
off rapidly away from midrapidity. It is suggested that 
this mid-rapidity peaking of the energy density may be 
related to the observed elliptical flow signal, which ex- 
hibit a similar pseudorapidity dependence. 
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XII. APPENDIX 



It is convenient to factorize the distribution 



(PN/drjdpT such that 



^<^^ dN,{r],pT) 

drjdpT dij dpT 



(14) 



Furthermore it is assumed that the relative abundance 
of pions, kaons, and protons/anti-protons is fixed as a 
function of pseudo-rapidity, i.e. 



dNi _ pdN 
drj * drj 

Under this assumption we may write 
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which is well supported by recent measurements of the 
identified particle distributions over a wide range of ra- 
pidity llUll. As illustrated in Fig. 7a, it is, however, 
clear that dependence on the rest mass of the particle 
is important only in the mid-rapidity region, which con- 
tributes only a relatively little to the total energy. 

The present analysis is based on the expression given 
in Eq. 6, which furthermore assumes that the average 
of dE/drj distributions is well approximated by a single 
distribution substituting the average particle mass, (m) 
obtained from the relative abundance of pions, kaons, 
and protons/ anti-protons. The fact that this is an ac- 
curate approximation is demonstrated in Fig. 7b, where 
the solid curve was obtained under this assumption. It 
is seen to coincide with the dashed curve representing a 
calculation using the expression of Eq. 5 with individual 
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values of the particle mass and mean transverse momen- 
tum of (pt)^=0.488 GeV/c, {pt)k=0-777 GeV/c, and 
(pt)p/p=1-126 GeV/c, which were obtained from con- 
verting the values [l8j measured in rapidity space to 
pseudo-rapidity. The success of this approximation is 
caused partly because of the relatively small abundance 
of kaons and protons, which is however known experimen- 
tally in the mid-rapidity region, where the total particle 
energy is dependent on the rest mass. At larger pseudo- 
rapidities, which account for the bulk of the energy, the 
results are insensitive to the particle masses. 



FIG. 8: Panel a): Comparison of dE/dr] calculated using a 
actual particle masses and the corresponding values of {pT)i 
(dashed curve) with one that uses an average particle mass 
(solid curve). A calculation with m = (dash-dotted curve) 
gives slightly smaller values near 77 = only. Panel b): Only 
in the midrapidity does the particle mass affect the energy; 
for I?;! > 3 the particle mass is irrelevant. 
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